QND measurements and state preparation in quantum gases by light detection 
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We consider light scattering from ultracold quantum gas in optical lattices into a cavity. The 
measurement of photons leaking out the cavity enables a quantum nondemolition (QND) access to 
various atomic variables. The time resolved light detection projects the motional state to various 
atom-number squeezed and macroscopic superposition states that strongly depend on the geometry. 
Modifications of the atomic and light properties at a single quantum trajectory are demonstrated. 
The quantum structure of final states can be revealed by ] 
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Quantum gases in optical lattices are of fundamental 
interest, as they provide an excellent testbed to study 
multipartite entanglement and many-body states, useful 
in condensed matter and quantum information [1]. Usu- 
ally, the role of light is reduced to a classical auxiliary 
tool for creating intriguing atomic states. In contrast, 
here we consider an ultimate quantum level, where quan- 
tum natures of both matter and light play a key role. 
This emerging level joining quantum optics [especially 
cavity quantum electrodynamics (QED)] and quantum 
gases, only recently became achievable and stimulated 
novel experimental Q and theoretical studies [1, 0, HI • 

We show that the atom-light entanglement enriches 
physics and enables the quantum nondemolition (QND) 
measurement and manipulation of atomic states. Ob- 
serving light allows to prepare different types of atom 
number squeezed and macroscopic superposition states. 
Note that the type of many-body states depends on the 
optical geometry. The quantum structure of final states 
can be revealed by further observations of the same sam- 
ple, which is an advantage over destructive schemes 

As we consider off-resonant interaction, independent 
of a particle level structure, our model might be also ap- 
plied to other phenomena in molecular physics [7] , where 
the molecule number fluctuations are important, and 
solid-state systems as semiconductors [Bose-Einstein con- 
densates (BEG) of exciton-polaritons] Q and supercon- 
ductors |y| (circuit cavity QED). Besides, the squeezed 
and macroscopic superposition states find applications 
in quantum interferometry and metrology [10]. 

Model. We consider (cf. Fig.[T]) N ultracold atoms in 
an optical lattice of M sites formed by strong off-resonant 
laser beams. A region of < M sites is also illuminated 
by a weak probe, which is scattered into a cavity. We will 
investigate, how the measurement of photons leaking the 
cavity will affect the atomic quantum state. 

The theory is based on the generalized Bose-Hubbard 
model taking into account the light quantization [1, Q . 
In contrast to Ref. we assume dynamics and mea- 
surement of the cavity mode is faster than atomic tun- 
nelins; Thus, neglecting the influence of tunneling 
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FIG. 1: (Color online) Setup. A lattice is illuminated by 
the transverse probe ao and probe through a mirror 77. The 
photodetector measures photons leaking the cavity. Due to 
the quantum back-action, the Ught measurement leads to the 
modification of the atomic quantum state. 



on light, we get the effective many-body Hamiltonian: 

H = Ti{ui + Uiibii)a[ai + fi/7io(^to<^o<^i + ^io<^o«i) 

—ih{r]*ai — 77aJ),(l) 

where ai is the cavity-mode annihilation operator and 
ao is the c-number probe amplitude of the frequencies 
LOi^p and spatial mode functions iii,o(r). Uim = gWm/^a 
(/, m = 0, 1), where ^1^0 are the atom-light coupling con- 
stants, Aa = coi—coa IS the cavity-atom detuning, r] is the 
probe through a mirror at uOp. We assumed the probe- 
cavity detuning Ap = uop — uoi <C A^. The operators 
Dim = Xlj=i ^i{^3)^rn{^j)T^j sum coutributious from all 
illuminated sites with the atom-number operators fij at 
the position Vj. Here, the rotating- wave approximation 
and adiabatic elimination of the exited state were used. 

The first term in Eq. ([T]) describes the atom-induced 
shift of the cavity resonance. The second one reflects 
scattering (diffraction) of the probe into a cavity. For a 
quantum gas the frequency shift and probe-cavity cou- 
pling coefficient are operators, which leads to different 
light scattering from various atomic quantum states [1] . 

The Hamiltonian ([T]) describes QND measurements of 
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variables related to Dim measuring the photon number 
a\ai fn\. Note, that one has a QND access to various 
many-body variables, as Dim strongly depend on the lat- 
tice and light geometry via iio,i(r). This is an advantage 
of the lattice comparing to single- or double-well setups, 
where the photon measurement back-action was consid- 
ered [l2[. Moreover, such a geometrical approach can be 
extended to other quantum arrays, e.g., ion strings p^ . 

For example, Du can reduce to the operator Nk of 
atom number at K sites [3]. If the probe and cavity 
modes are coupled at a diffraction maximum (Bragg an- 
gle), i.e., all atoms scatter light in phase, the probe-cavity 
coupling is maximized, Dio = Nk- If they are coupled 
at a diffraction minimum, i.e., neighboring atoms scatter 
out of phase, Dio = Y^f=i{—'^y~^^nj is the operator of 
number difference between odd and even sites. Thus, the 
atom number as well as number difference can be non- 
destructively measured. Note, that those are just two 
of many examples of how a QND- variable, and thus the 
projected state, can be chosen by the geometry. 

Measurement back-action. The expression for the 
initial motional state of atoms reads 



l*(0)) 



q 



(2) 



which is a superposition of Fock states reflecting all possi- 
ble classical configurations ^ = {<7i, .., Qm} of N atoms at 
M sites, where qj is the atom number at the site j. This 
superposition displays the uncertainty principle, stating 
that even a single atom can be delocalized in space. 

While interacting, the light and atoms get entan- 
gled. Quantum mechanics predicts that measurements of 
one subsystems (light) provides conditional information 
about, or affects, another one (gas). We will show, how 
the atomic uncertainty is affected by the light detection. 

We use the open system approach Q for counting 
photons leaking the cavity of decay rate hi. When a pho- 
ton is detected, the jump operator is applied to the state: 
|^c(^)) <^i|^c(^))- Between the counts, the system 
evolves with a non-Hermitian Hamiltonian H — ihtva\ai. 
Such an evolution gives a quantum trajectory for |^c(^)) 
conditioned on the detection of photons at times ti, ^2, ... 

It is known [15] that, if a coherent probe illuminates a 
classical atomic configuration in a cavity, the light state is 
proportional to a coherent state \aq{t)) with aq{t) given 
by the classical Maxwell's equation. Thanks to the ap- 
proximation, where the tunneling does not affect light, we 
can get a simple analytical solution of the coupled light- 
matter dynamics. Each atomic Fock state in Eq. ([2|) will 
be correlated with a coherent light state with param- 
eters given only by the corresponding configuration q: 
l^c(t)) = EqC^q^MMt)]\QU".QM)\aq{t))/F{t), whcrc 
F(t) gives the normalization. So, the problem to find 
|^c(^)) reduces to finding aq{t), ^q(^) for all classical 
configurations forming the initial |^(0)). Although a so- 
lution is available for any we present it for t > 



when the steady state is achieved in all ag(t), and as- 
suming the first photon was detected at ti > 

Due to the steady state in all ^^(t), the solution is 
independent of the detection times and after m counts is 

\^c{m,t)) = J-S2a^e^^(''>4\qu...,qM)K), (3) 

^ T] - iUioaoDl^ 

z([/nI^?i-A^) + ^' ^ ^ 
^q{t) = -\aq\^i^t + {r]al - iUiQaoDl^al - c.c.)t/2, (5) 

where D^^ = Yl/f=i^i{^3)^ra{^j)(lj is a realization of 
Dim at {^1, ..^m}; <^o, ^5 and aq all oscillating in steady 
state at uOp were replaced by their constant amplitudes. 

As we see, each light amplitude ocqit)^ Eq. (jH, is given 
by a Lorentzian corresponding to classical optics. Eq. (|3|) 
shows that the probability to find an atom configuration 
q, pq{m,t) = l^gp^ exp (— 2|agp/^t)|c^p/F^, changes in 
time due to the photodetection. This demonstrates the 
back-action of the light measurement on the atomic state. 

In the following, we will show consequences of Eq. ([3]) 
for two cases, where only one probe (ao or 77) exists. For 
transverse probing (ao 7^ 0), we also neglect the mode 
shift, assuming UuDf-^ <C hz or Ap. Thus, in both ex- 
amples, aq dl]) depends on the configuration q only via a 
single statistical quantity now called z: z = D\^ for cav- 
ity probing (77 7^ 0), and z = D\q for transverse probing. 

From Eq. ([3|) we can determine the probability distri- 
bution of finding a given z after time t as 



p{z,m,t) 



\2m -2\az\'^Kt 



(6) 



where the initial distribution po{z) = Xlg' kg'P^ such 
that all configurations q^ have the same = 
\^z\'^^ ^'^V {—'^\<^z\'^ f^t)Po{z) provides normalization. 
Transverse probing at diffraction maximum. As 
was mentioned, at the Bragg angle, Dio = Nk is the 
operator of atom number at K sites. So, z varies from 
to N reflecting possibilities to find any atom number at K 
sites. The light amplitudes ^ az = Cz are proportional 
to C = iUioao/{iAp — hz). The probability (|6]) reads 



p{z,m,t) = z''^e-' ^Po{z)/F' 



(7) 



with a characteristic time r = 2|Cp/^t. 

When time progresses, both m and r increase with 
a probabilistic relation between them. The Quantum 
Monte Carlo method [13] establishes such a relation, thus 
giving a trajectory. Note, that thanks to the simple an- 
alytical solution (|3|), it gets extremely simple. In each 
step, it consists in the calculation of the photon num- 
ber in the state given by Eq. (|3]) and comparing it with 
a random number generated in advance, thus, deciding 
whether the detection or no-count process has happened. 

If the initial atom number z at K sites is uncer- 
tain, po{z) is broad [for the superfluid (SF) it is nearly 
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FIG. 2: (Color online) Photodetections at diffraction max- 
imum, (a) Shrinking atom number distribution at different 
times r = 0, 0.005, 0.018, 0.03, 0.05, 0.5 (A-F); (b) decreas- 
ing width (c) stabilizing mean atom number {z)c' Initial 
state: SF, N = 100 atoms, K = M/2 = 50 illuminated sites. 

Gaussian [3]], and Eq. ([7j) shows that p{z^m^t) is 
strongly modified during the measurement. The function 
z'^^ exp {—z^r) has its maximum at Zi = \/m/r and full 
width at half maximum (FWHM) 5z ^ V21n2/r (for 
Sz <C zi). Thus, multiplying po{z) by this function will 
shrink the distribution p(z, m, t) to a narrow peak at Zi 
with the width decreasing in time (Fig. [2j). 

This describes the projection of the atomic quantum 
state to a final state with squeezed atom number at K 
sites (a Fock states \zi, N — zi) with zi atoms at K sites 
and N — zi atoms at M — i^ sites). When Sz < 1^ the final 
collapse is even faster than y^, due to the discreteness 
of p{z^ m, t). Measuring the photon number m and time 
t, one can determine zi of a quantum trajectory. 

In contrast to recent results in spin squeezing [l6t. 
which can be also obtained for thermal atoms [171], in 
our work, quantum nature of ultracold atoms is crucial, 
as we deal with the atom number fluctuations appearing 
due to the delocalization of ultracold atoms in space. 

After the distribution shrinks to a single Zi, the light 
collapses to a single coherent state l^^i), and the atoms 
and light get disentangled with a factorized state 

|*c) = |^i,iV-zi)K). (8) 

So, light statistics evolves from super-Poissonian to Pois- 
sonian. The conditioned (i.e., at a single trajectory) cav- 
ity photon number {a[ai)c{t) = |Cp ^^=o z'^p{z^m,t) is 
given by the second moment of p{z,m,t). Its dynamics 
[very similar to {z)c in Fig.[2l[c)] has jumps, even though 
all az{t) are continuous. In the no-count process, {a\ai)c 
decreases, while at one-count it jumps upwards, which is 
a signature of super-Poissonian statistics. Finally, it re- 
duces to {a\ai)c = reflecting a direct correspon- 
dence between the final atom number and cavity photon 
number, which is useful for experiments. 
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FIG. 3: (Color online) Photodetections at diffraction mini- 
mum, (a) Shrinking distribution of the atom-number differ- 
ence for various times r = 0, 0.02, 0.03, 0.17, 0.65 (A-E). The 
doublet corresponds to macroscopic superposition state, (b) 
Relative photon number {a\ai)c/\C\'^ with quantum jumps. 
Initial state: SF, N = 100 atoms, K = M = 100 sites. 



Even the final Fock state still contains the atom-atom 
entanglement, as many components |(7i,..,^m) can have 
the same zi. For example, the SF state can be rep- 
resented as \SF)n,m = T.z\^\SF)z^k\SF)n-z,m-k 
{Bz are binomial coefficients). After the measurement, it 
ends up in \SF) z:^,k\SF) n-zi,m-k ^ i.e., the product of 
two uncorrelated superfluids. 

Our measurement scheme determines (by squeezing) 
the atom number at a particular lattice region and 
projects the initial atomic state to some subspace. How- 
ever, the atom number at different regions keeps quantum 
uncertainty. So, the quantum structure of the final state 
can be revealed in a further optical or matter- wave exper- 
iment. Thanks to the lattice geometry, one can change 
the illuminated region, and further study measurement- 
induced collapse of the state in the remaining subspace. 

Even in matter- wave experiments [l|, the product of 
SFs will look different from the initial SF: the atoms from 
different regions will not interfere in average. Note, that 
we did not specify how K sites were selected. One can 
illuminate a continuous region. However, one can illumi- 
nate each second site by choosing the probe wavelength 
twice as lattice period and get number squeezing at odd 
and even sites. In this way, one gets a measurement- 
prepared product of two SFs "loaded" at sites one by 
one (e.g. atoms at odd sites belong to one SF, while at 
even sites to another). While the initial SF shows the 
long-range coherence (blbj) with the lattice period, the 
prepared state will demonstrate the doubled period in 
(blbj) {bj is the atom annihilation operator). 

Transverse probing at diffraction minimum. In 
contrast to classical atoms, quantum gases scatter light 
even in diffraction minima [3]. Here z = D\q = 
^j=i{—^y~^^Qj is the atom number difference between 
odd and even sites, varying from —N to N with a step 
2 (we assumed K = M). Eq. ([7j) keeps its form with a 
new meaning of z and Pq{z) [for SF, new po(^) is nearly 
a Gaussian centered at 2: = and the width [3]]. 

The striking difference from the diffraction maximum 
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is that our measurement ([7j) is not sensitive to the sign of 
while the amphtudes = Cz are. So, the final state 
is a macroscopic superpositions of two Fock states with 
Zi^2 = ±\An/V and light amplitudes: = — <^2i, 

|*c) = + (-in -z^)\- a^J)/V2. (9) 

Figure [3] shows the collapse to a doublet probability 
p(±zi,m,t) and the photon- number trajectory, where 
upward jumps and no-count decreases can be seen. 

In contrast to a maximum, even in the final state, the 
light and matter are not disentangled. Moreover, to keep 
the purity of the state, one should know precisely the 
number of detected photons, because of the sign fiip in 
Eq. ([9|). This refiects the fragility of macroscopic super- 
position states with respect to the decoherence. 

As a result, the measurement-based state preparation 
at the diffraction maximum ([8]) is much more robust. 

For probing through a mirror (77 7^ 0, ao = 0), in con- 
trast to the transverse probing, the probability distribu- 
tion can collapse both to a singlet and doublet. Impor- 
tantly, here the superposition state can be more robust 
than (|9]). This is due to a smaller phase jump between 
two Fock states, which can be obtained from Eqs. ([3])- ([5]). 

Deviations from our ideal setup can slightly modify the 
results, i.e., the mode profiles will lead, instead of atom 
number, to more general variables given by Dim- 

In summary, we showed that using the light-matter 
entanglement in ultracold gases enables QND measure- 
ments of different atomic variables and creation of spe- 
cific atomic states. The state type is determined by the 
optical geometry. Our model can be generalized to other 
quantum arrays. Cavity QED with quantum gases can 
operate with atom numbers from millions to one [13] • 
Thanks to recent experimental breakthroughs [2] , prepar- 
ing various kinds of atom number squeezing is already 
doable, and creation of the superposition states with, at 
least, small particle number [10| may become practical. 
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